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Abstract 

Phase diagrams of the two-dimensional one-band t-t' Hubbard model are 
obtained within the two-patch and the temperature-cutoff many-patch renor- 
malization group approach. At small t' and at van Hove band fillings anti- 
ferromagnetism dominates, while with increasing t' or changing filling anti- 
ferromagnetism is replaced by d-wave superconductivity. Near t' = t/2 and 
close to van Hove band fillings the system is unstable towards ferromagnetism. 
Away from van Hove band fillings this ferromagnetic instability is replaced by 
a region with dominating triplet p-wave superconducting correlations. The re- 
sults of the renormalization-group approach are compared with the mean-field 

results and the results of the T-matrix approximation. 
PACS numbers: 71.10.Fd; 71.27.-ha; 74.25. Dw 
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I. INTRODUCTION 



The close relation between antiferromagnetism (AF) and d-wecve superconductivity (dSC) 
was the subject of intensive investigations during the last two decades (see e.g. Refs. [1-5]). 
In particular, it was argued that the superconducting properties of high- Tc (HTSC) materi- 
als are intimately related to their inherent antiferromagnetic correlations and many features 
of these materials were explained from the point of view of competition between antiferro- 
magnetic and superconducting correlations [3] . On the other hand, AF spin fluctuations also 
serve as the natural candidate for the pairing mechanism of dSC [4,5]. A distinctly different 
physical situation is realized in the layered ruthenate Sr2Ru04, which is an unconvential 
and most likely triplet superconductor [6] . It was proposed, that the pairing in this material 
results from ferromagnetic spin fluctuations [7,8]. Although inelastic neutron scattering has 
so far been unsuccesful to detect significant low-energy ferromagnetic spin fluctuations in 
this material [9] , this idea finds experimental support from the recent measurements of the 
susceptibility of the electron doped compound Sr2-a;Laj;Ru04 [10] which revealed the ten- 
dency towards ferromagnetism with La doping. Furthermore, the isoelectronic compound 
Ca2Ru04 also shows ferromagnetism under hydrostatic pressure [11]. 

Both, copper-oxide systems and Sr2Ru04, are layered materials. Therefore both systems 
motivate the investigation of the competition and the mutual interplay between magnetic 
and superconducting instabilities in two-dimensional (2D) correlated electron systems. For 
this type of analysis it is important to account for specific band structure related phenomena, 
namely for the form of the Fermi surface (FS) and the electronic dispersion. The influence 
of the shape of the FS on superconducting and magnetic properties is of interest from both, 
theoretical and experimental points of view and the theoretical analysis can be guided by 
material-specific information obtained from angle-resolved photoemission (ARPES) experi- 
ments [12-15]. 

The simplest theoretical model which allows to investigate the effect of the band dis- 
persion on magnetic ordering and superconductivity of 2D systems is the single-band t-t' 
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Hubbard model on a square lattice which takes into account both nearest-neighbor t and 
next-nearest-neighbor t' hopping. This model is often discussed in connection with HTSC 
compounds, and it describes well the shape of the FSs of cuprate superconductors observed 
in ARPES [12-14]. In particular, the value t'/t — 0.15 was chosen for La2Cu04 and the value 
t'/t — 0.30 for the Bi2212 system [16] in the tight- binding parametrization of the relevant 
electronic band for the Cu02-planes, although the realistic modelling of the latter bilayer 
material requires the inclusion of interlayer hopping as well. On the other hand, Sr2Ru04 
has three relevant bands [17]. Interband effects are not negligible in this material, and may 
even prove important for the origin of unconventional superconductivity [18]. 

Already in early mean-field and quantum Monte Carlo (QMC) studies of the t-t' Hubbard 
model [19] it was found that depending on the ratio t'/t and the band fiUing, different types 
of instabilities are possible. For small t'/t near half-filling the FS is almost nested, which 
is the origin for antiferromagnetism in the weak-coupling regime, t' hopping destroys the 
perfect nesting property of the FS and therefore leads to "frustration" of antiferromagnetism 
due to the hopping processes on the same sublattice and may therefore favor the emergence 
of a superconducting state [20] . Furthermore, t' hopping also weakens the tendency towards 
stripe formation [21] and by the suppression of this alternative instability superconducting 
fluctuations may get enhanced. At the same time, larger values of t' move the system closer 
to a ferromagnetic instability, since for t'/t close to 1/2 the dispersion is flattened close to 
the bottom of the band. This leads to fiat-band ferromagnetism [22] at low densities which 
was investigated earlier for the 2D t-t' Hubbard model within the T-matrix approximation 
[23,24] and projected QMC simulations [25]. 

The interplay of antiferromagnetism and d-wecve superconductivity in the one-band t-t' 
Hubbard model was recently reconsidered within many-patch renormalization-group (RG) 
approaches [26-28]. It was indeed verified that with increasing t' and/or decreasing band fill- 
ing antiferromagnetism is replaced by d-wave superconductivity. In the early RG approaches 
of Refs. [26-28] particle-hole scattering at small momenta was not treated on equal footing 
with other types of scattering, and therefore these analyses were unable to search for a pos- 
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sible ferromagnetic instability (see the discussion in Refs. [29,30]). It was shown in Ref. [31] 
however, that particle-hole scattering at small momenta does indeed lead to the appearance 
of a ferromagnetic phase at large enough t' /t and at van Hove (vH) band fillings; the onset of 
dominant ferromagnetic correlations was found to occur for t' /t > 0.27. However, unlike in 
Refs. [26-28], the contribution of the Cooper channel was not taken into account in Rcf. [31]. 
The possibility of a ferromagnetic instability was also investigated within a simplified two- 
patch RG scheme [29] , which considers only the scattering of electrons in the vicinity of the 
"singular" points (tt, 0) and (0, n) in momentum space and therefore gives only a rough pic- 
ture for the RG scaling behavior of the coupling constants. The temperature-cutoff version 
of the many-patch RG approach (TCRG) recently introduced by Honerkamp and Salmhofer 
[30] includes the contributions of the whole Brillouin zone and uses the temperature as a 
natural low-energy cutoff parameter in order to avoid the technical difficulties with the inclu- 
sion of small-momentum particle-hole scattering. It was demonstrated in Refs. [29,30] that 
the proper account of all scattering channels indeed leads to ferromagnetism at large enough 
t'/t. Moreover, the critical value of it' /t)c for the stability of ferromagnetism is ^/-dependent 
[29] (unlike the results of Ref. [31]), in particular for [/ — > ferromagnetism exists only in 
the flat-band low-density limit {t'/t)c — > 1/2 in qualitative agreement with the results of 
the T- matrix approximation for the effective electron-electron interaction vertex [23]. Nat- 
urally, ferromagnetic and d-wave superconducting fluctuations tend to suppress each other 
[29,30]. The suppression of ferromagnetism by superconducting fluctuations is reminiscent 
of the well-known Kanamori screening [32] . On the other hand, as shown in Ref. [30] , the 
tendency towards triplet superconductivity is enhanced by ferromagnetic fluctuations and 
may exist in the vicinity of a ferromagnetic phase. Note that the Pomeranchuk instability, 
which was proposed for small t' in Ref. [33], was shown to be a non-leading instability in 
the t-t' Hubbard model [34]. 

With these recent results it appears as a natural task to investigate systematically the 
weak-coupling phase diagram of the 2D Hubbard model within the RG approach. Earlier, 
such an analysis was performed only with the momentum cutoff RG versions [26-28] , which, 
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as we discussed above, do not allow to include the contribution of ferromagnetic fluctuations. 
Previous studies of the model with the inclusion of all electron scattering channels were 
performed either in the two-patch RG scheme [29], which is restricted to vH band fillings 
or within the TCRG approach [30,34], which however was apphed only for some selected 
parameter values and therefore results for a phase diagram in the weak-coupling regime were 
not obtained. 

The purpose of the present paper is to investigate systematically different types of insta- 
bilities of the t-t' Hubbard model within two-patch and the many-patch TCRG approach. 
The paper is organized as follows. In Section II we give a summary of the RG methods we 
use. In Section III we present the phase diagrams obtained and compare the results with 
previous investigations of the t — t' Hubbard model. In Section IV we discuss the results 
and conclude. 

II. RENORMALIZATION-GROUP APPROACHES 

We consider the t-t' Hubbard model 

H^-Yl iiAa(^j- + ^ 51! "^T^a - - (1) 

ija i 

where the hopping amplitude tij — t for nearest neighbor sites i and j and tij — —t' for next- 
nearest neighbor sites (t, t' > 0) on a square lattice (we have shifted the chemical potential 
by At' for further convenience). In momentum space Eq. (1) reads 

H = ^e-kci^C],„ + ^ ^ Xl^ki<T42(T'Ck3a'Ck4a5ki+k2-k3-k4 (2) 
kcr kik2k3k4 era' 

where the Kronecker 5-symbol ensures momentum conservation and the dispersion has the 
form 

£k = — 2t(cos kx + cos ky) -\- 4t'(cos k^ cos ky + 1) — fi (3) 

where the lattice constant is set to unity. The tight-binding spectrum (3) leads to vH 
singularities (vHS) in the density of states arising from the contributions around the points 
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— (tt, 0) and = (0,7r). These singularities lie at the FS if = 0. For t' — the 
corresponding filling is uvh — 1 and the FS is nested, but the nesting is removed for any 
t' ^ 0. The dependence of the vH band filling on t' is shown in Fig.l, and the shape of the 
FS at different t'/t and vH band fillings is shown in Fig. 2. 

The standard RG strategy for fermion systems [35] is to integrate out step by step the 
electronic states which are far from the FS (i.e. the states with the energy A — dA < 
£k < A at each RG step). This procedure meets a difficulty when it is applied to a FS 
with singular points, i.e. the points kp with vanishing Fermi velocity V£:k|k=k|, = 0, as 
in points k^ and kg. In this case, the states with the same excitation energy £k become 
inequivalent: the excitations with momenta closer to the singular points produce more 
divergent contributions to the renormalization of the electron-electron interaction vertices 
than the excitations with momenta far from the singular points. Therefore, an additional 
separation of the momenta besides the standard separation into "slow" {e^ < A) and "fast" 
(sk > A) modes is needed. The two-patch approach which we consider in Sect. Ilia accounts 
only for the most singular contributions coming from the immediate vicinities of the singular 
points. The more sophisticated many-patch approaches of Refs. [26-28,30] (see Sect. Illb) 
take into account the momentum dependence of the interaction in a more accurate way 
by introducing a set of patches which cover the entire Brillouin zone and parametrize the 
interactions by the position of incoming and outgoing momenta on the patched FS. 

A. Two-patch renormalization-group approach 

The two-patch approach [36,37,29] is restricted to the vH band fillings only. At these 
fillings the density of states at the Fermi energy and the electron-electron interaction vertices 
at momenta k, = k^ ^ contain logarithmic divergencies arising from the momentum integra- 
tions in the vicinity of the points k^^^. Therefore these contributions are the most important 
for the calculation of the renormalized interaction vertices. Accordingly, we subdivide the 
momentum space into three types of regions (see Fig. 3). Regions I with k e 0(^4) V 0{B) 
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where 

0{A) = {k : |k - k^l < A A \ek/t\ > e-^/\^-^^\} (4) 

and similar for 0{B) (A is a momentum cutoff parameter) produce the most singular contri- 
bution to the renormalization of the vertices. Regions II contain the electronic states which 
are close to the FS but far from vH singularities. It can be proven that the contributions 
of regions II to the renormalization of the vertices is subleading in comparison with the 
contributions of regions I, provided that t'/t is not small, i.e. if the nesting effects are not 
important. Finally, regions III contain the excitations which are far from both, the FS and 
vHS and do not produce diverging contributions to any quantity. Therefore, in the simplest 
approximation it is reasonable to neglect the contributions of regions II and III altogether. 
A more accurate treatment within the many-patch RG approach will be performed in the 
next section. 

To account for the excitations with momenta in regions I, it is convenient to introduce 
new electron operators Ok and 6k by 

I ak-k^,a k e 0(A) 

Cka = < 

[ /^k-ks,. k e 0{B) 

For momenta k e 0(^4) V 0{B) in the vicinity of the vH points the dispersion is expanded 
as 

£k^+p = £p = -2t(sin^ V^pI- cos^ - (5a) 
= -Itp+p- - II 

£kB+p = ^p = 2i(cos^ (ppl - sin^ (ppD - (5b) 
= 2tp+p^ - /X 

where cos{2(fi) — R — 2t' / 1, p± — Px sin (fi^Py cos (p, and p± — Px cos (p^Py sin (p. Using the 
new electron operators we write the effective Hamiltonian in the form 

pa pa 
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Pi, era' 

(6) 

Pi,tTa' 

where 

A = ln(A/ max(pj+,pi_,pj+,pi_, T/t); (7) 

the summation in Eq. (6) is restricted to momenta Pi with |pi| < A and I^^a B+Pi/A > 
e-A/|Pil. 

As shown in Fig. 4, the vertices gi to 514 represent different types of scattering processes 
of electrons with momenta close to the vHS. The bare value for all four vertices is = 
U I [A.'kH^ . The momentum dependence of the vertex inside regions I is accounted for through 
the scaling variable A only. Note however, that the momentum dependence of the electronic 
spectrum within each patch is correctly taken into account in the two-patch approach. 

To obtain the dependence of the vertices gi on A we integrate out at each RG step 
the fermions a-p with momenta Ae~^ < P± < Ae~^~''^, and fermions 6p with momenta 
Ae"'*' < p± < Ae~'^~'^'^ (see the detailed description in Ref. [38]). As argued above we 
neglect the renormalization of the gi arising from the regions II and III of Fig. 3 since this 
leads to subleading corrections at weak coupling. We determine the RG equations for the 
vertices gi{X) in the form [36-38,29] 

dg'i/dA = 2di{X)gi{g2 - gi) + 2d2gm - 2 d3gig2, 
dg2/dX = di{X){gl + gl) + 2^2 (5-1 - g2)94 - d^i^gl + gl), 
dg'a/dA = -2do{X)g3g4 + 2di{X)g3{2g2 - gi), 

dg^/dX = -do{X){gl + gl) + 0^2 (^1' + 2^i^2 - 2^2' + qI), (8) 

where 

do{X) = 2A/V1 - d2 = 2/Vl-R^; 

= 2tan-\R/Vl- R'^)/R; 
cii(A) = 2min{A,ln[(l + vT^rR2)/i?]}. (9) 



Eqs. (8) have to be solved with the initial conditions gi{0) — gf. 

The authors of Ref. [31] argued that the kinematic restrictions lead to the absence of 
particle-particle scattering contributions to the vertices gi {do = = in our notations). 
This conclusion however is connected with the difficulty of the infinitesimal version of Wil- 
son's RG approach with a sharp momentum cutoff, since this approach does not allow to 
treat correctly the vertices with nonzero momentum transfer and gives artificially no renor- 
malization for such vertices [39] . That is why in the present approach we consider the vertices 
at the special vH points only, rather than considering vertices with arbitrary momentum 
transfers. Note that this difficulty does not arise in Wilson's RG approach with a smooth 
cutoff and/or discrete RG transformations [39]; in this case particle-particle scattering does 
contribute to the renormalization of the vertices with arbitrary momenta. 

In order to explore the possible instabilities of the system, we consider the behavior of 
the zero-frequency, time-ordered response functions 

l/T 

Xm^ J dT{T[dliT)d^m) (10) 


in the zero-temperature limit T ^ 0. Om(r) are the following operators 

k 

OdSC =j;^Yl /k(^4,aClk,-a' 
kfT 

= -^Xl^^U^k,^, (11) 
ko- 

in the Heisenberg representation, T is the imaginary time ordering operation, and Q = 
(tt, tt). The order parameters which correspond to p-wave pairing 

^ps^c^^E^^'^i^^^V.' (12) 

ko- 

with h^^ — sin kx,y are irrelevant with the restriction of momenta to the vicinities of vH 
points, since /ik^ b — ^ therefore the possibility of triplet pairing can not be considered 
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within the two-patch approach. However, these order parameters can be taken into account 
in many-patch approaches (see Section Illb). 

Picking up the logarithmical divergences in A we obtain the RG equations for the dimen- 
sionless susceptibihties Xm — '^^'^'tXm in the same approximations as discussed above (cf. 
Refs. [36,37,29]): 

dx„(A)/dA = 4^(A)7^^(A), (13) 
dln7e„(A)/dA = 4^(A)r^(A), 

where the coefficients (m — AF, dSC, or F) are given by 

Taf = 92 + 93] Tp = gi + g^; 

(14) 

Tdsc = 93~ 94- 

In Eqs. (13) a^sc = 0, caf = 1; and ap = 2. Eqs. (13) have to be solved with the initial 
conditions TZmi^) = 1, Xm(0) = 0. 

The numerical solutions of Eqs. (8) show, that at a critical value Ac of the scahng 
parameter A some of the vertices and susceptibihties are divergent. For a given Ac the size A 
of the patches is restricted by ln(4/A) <^ Ac. The latter criterion follows from the condition 
that the contribution of the electrons with \k±\ < A to particle-hole and particle-particle 
bubbles is dominant (sec e.g. Ref. [40]). Wc choose A = 1 and require Ac ^ In 4 ~ 1. Since 
Ac decreases with increasing interaction strength, this criterion defines the interaction range 
where the two-patch RG approach is valid. 

As an example, we show in Figs. 5,6 the result of the numerical solutions of Eqs. (8) for 
U — 2t and two different choises of t'/t — 0.1 and t' /t — 0.45. The behavior of the coupling 
constants is qualitatively different in these two cases. While in the first case we have g2,3 
flowing to +00, gi to — oo, and gi is mostly unchanged during the RG flow (we denote the 
corresponding combination {rn-\- -\ — ), the signs correspond to the behavior of the coupling 
constants g^-g^, m means marginal. In the second case 511,2,4 grow to -|-oo while g^ goes to 
zero, i.e. we observe (-1- + 0+) behavior of the coupling constants. The comparison of the 
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corresponding susceptibilities shows that for t' /t = 0.1 the antiferromagnetic susceptibility is 
the most divergent, while in the case t' jt = 0.45 the ferromagnetic susceptibility dominates 
and therefore the two different coupling constant flows reflect two different instabilities of 
the system. We discuss the complete phase diagram at vH band fillings in Sect. Ilia. 

B. Many-patch renormalization-group analysis 

In the many-patch analysis we follow the temperature-cutoff RG for one-particle irre- 
ducible Green functions proposed recently by Honerkamp and Salmhofer in Ref. [30]. This 
version of the RG uses the temperature as a natural cutoff parameter, allowing to account 
for both the excitations with momenta close to the FS and far from it, which is necessary 
for the description of instabilities which arise from zero-momentum particle-hole scatter- 
ing, e.g. ferromagnetism. Neglecting the frequency dependence of the vertices, which is 
considered not important in the weak-coupling regime, the RG differential equation for the 
temperature- and momentum-dependent electron-electron interaction vertex has the form 
[30] (see Fig. 7) 

^\/T(ki, k2, kg) = - ;^ XI V^T(ki, ka, p)Lpp(p, -p + ki + k2)yT(p, -p + ki + ka, kg) 

p 

+^ J] [ ~ 2yr(ki, p, k3)yT(p + ki - kg, ka, p) + ^rlki, p, k3)yr(k2, p + ki - kg, p) 
p 

-M^T(ki, p, p ki - kg) Vt(p + ki - kg, k2, p)]Lph(p, p + ki - kg) 

+ X yr(ki, p + k2 - kg, p)Lph(p, p + ks - k3)yT(p, ks, kg) (15) 
p 

where 

/^(£k) - /^(£k') 



^ph(k,k') 



Lpp(k,k') = :^i^^^i±^^, (16) 
£k + £k' 

and /^(e) = d/(£)/dT', f{e) is the Fermi function. Eq. (15) has to be solved with the 
initial condition VTo(ki, k2, kg) = U where the initial temperature Tq is of the order of the 
bandwidth. 
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The evolution of the vertices with decreasing temperature determines the temperature 
dependence of the susceptibihties according to [30,34] 

-^XmT = ^7^„T(k')7^„^T(Tk' + q^)Lpp,ph(k', ^k' + q„), (17) 
k' 

-^7^^^(k) ^mT(k')r^r(k, k')Lpp,ph(k', ifk' + q^) 



dT 

k' 



where 



yT(k,k',k' + qJ form = AForF 
^T(k,k) = . (18) 

VriK -k, k') for m = dSC or pSC 

Qm — Q for the AF susceptibihty and q^ = otherwise. The upper signs and pp indices in 
Eq. (17) refer to the superconducting instabihties (dSC and pSC), and the lower signs and 
ph indices to the other susceptibilities. The initial conditions for Eqs. (17) are 



cos kx — cos ky for dSC 



7em,To(k) = { 



sin kx,y for pSC , (19) 
1 otherwise 



and Xm,To — 0. To solve numerically Eqs. (15) and (17), we use the discretization of momen- 
tum space in Np — 48 patches and the same patching scheme as proposed in Ref. [30]. By 
exploiting the symmetries of the square lattice, this reduces the above integro-differential 
equations (15) and (17) to a set of 5824 differential equations which were solved numeri- 
cally. We use the value of the starting temperature To = 12i, which is slightly larger than 
the bandwidth, and we stop the flow of the coupling constants at the temperature Tx when 
the maximum absolute value of the vertex function is larger than \4iax — 18t. Note that 
the initial k-dependence of the vertices (19) is slightly changed during the RG flow: the 
vertices which correspond to d- and p-wave superconductivity acquire g- and /-wave (and 
even higher order) harmonics, respectively, and the vertices which have s-wave symmetry 
(F, AF) acquire an additional extended s-wave {cos k^ + cos ky) component. However, these 
additional corrections are small. 
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To extract more detailed information about different instabilities we extrapolate the 
inverse susceptibilities XmT temperatures lower than Tx- For magnetically ordered or 
superconducting ground states the corresponding extrapolated inverse susceptibilities 
vanish at some temperature below Tx- The vanishing of the inverse susceptibilities at 
finite temperatures is an artifact of the one-loop RG approach and should be understood to 
determine a crossover temperature into a renormalized classical regime with strong magnetic 
or superconducting fluctuations and exponentially large correlation length ^ oc ex.p{A/T) 
for F and AF instabilities (see Appendix) and $, oc exp{A/^/T'^^^BKT) for superconduct- 
ing instabilities. Tbkt is the Berezinskii-Kosterlitz-Thouless transition temperature which 
arises because the superconducting transition in 2D is in the same symmetry class as the 
classical 2D XY model [44,45]. Although the inverse magnetic susceptibilities must be finite 
at T < Tp^^p, they are exponentially small, since in this regime Xf af R^fs. 
[41-43]. The same concerns the superconducting susceptibilities in the temperature range 
Tbkt < T < T^sc.pSC' where critical behavior Xdsc psc °^ (77 = 1/4 is the critical 

exponent for the susceptibility) as for the XY model [44] is expected. At T < Tbkt the su- 
perconducting correlations have power law decay in real space and the inverse static uniform 
order parameter susceptibility does indeed vanish. 

III. RESULTS AND PHASE DIAGRAMS 

A. van Hove band fillings 

First we trace the vH band fillings at different t', which are determined by the condition 
IJ, = 0. The corresponding phase diagram in t' — U coordinates is plotted in Fig. 8. Solid 
lines correspond to the phase boundaries obtained within the two-patch approach, while the 
symbols show different types of instabilities obtained within the many-patch RG scheme 
[46]. Our many-patch results for U — 3t agree quantatively with those obtained previously 
in Ref. [30]. 
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In both approaches we find antiferromagnetism for small t' , ferromagnetism for t' close to 
1/2 and c?-wave superconductivity for intermediate t' . The t' range with a tendency towards 
(iSC decreases with increasing U. For intermediate t' /t the susceptibilities with momenta 
Q =(7r — 5, tt) are stronger than the antiferromagnetic susceptibility. The incommensurate 
magnetic regions are indicated in Fig. 8 as well. In the region t' /t ~ 0.3 the behavior of the 
coupling constants in the two-patch approach becomes "frustrated" , namely all gi ^ with 
increasing A. This frustration is the consequence of the competition of antiferromagnetic 
and superconducting instabilities from one side and the ferromagnetic instability from the 
other side and therefore at t' /t > 0.2 neither the antiferromagnetic nor the superconducting 
susceptibility diverges in the two-patch approach. The many-patch approach suffers less 
from this problem; the frustrated behavior of the vertices is observed only very close to the 
boundary of ferromagnetic and antiferromagnetic or superconducting phases. 

The boundary to the ferromagnetic phase appears almost identical in two- and many- 
patch approaches, but the two-patch approach fails to reproduce the location of the phase 
boundary between the antiferromagnetic and superconducting phases for small t'. This is 
similar to the results for the extended U-V-J Hubbard model [38]. As mentioned in Section 
Ila, this difference is traced to the same behavior (m + -\ — ) of the coupling constants 
on approaching AF and dSC instabilities, while the ferromagnetic phase is signalled by a 
different behavior of the coupling constants (-I--I-0-I-). Furthermore, the near-nesting effects, 
which are not accounted for in the two-patch RG approach become particularly important 
at small t'/t. 

We also mark in Fig. 8 the result of Alvarez et al. [31] for the boundary of the ferro- 
magnetic phase, obtained by neglecting the contribution of particle-particle scattering. In 
this case the corresponding two-patch RG equations can be solved analytically, since the 
coefficients di in Eq. (8) become A-independent. The main difference in comparison with 
Ref. [31] is that particle-particle scattering leads to a [/-dependence of the critical value 
{t'/t)c for the appearance of ferromagnetism so that {t'/t)c — > 1/2 for f/ ^ in qualitative 
agreement with the results of the T- matrix approximation [23]. At the same time, neglect- 
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ing particle-particle scattering gives {t'/t)c ^ 0.27 independent of U. Note that the value 
(t'/t)c is determined by the condition of equal non-interacting particle-hole susceptibilities 
Xo(0) = Xo(Q) and therefore coincides with the "mean-field" criterion for the boundary of 
the ferromagnetic phase. 

B. Antiferromagnetic instability at small t' 

The case t' = for fillings close to n = 1 was intensively studied previously within 
momentum-cutoff RG approaches [26,27]. We plot the results of TCRG a,t t' — in Fig. 
9. There is a hne of critical concentrations ndU) such that ior n < ndU) the extrapolated 
inverse susceptibility x^f (^oes not reach zero for any temperature and the ground state is 
expected not to have long-range antiferromagnetic order. For comparison, we also plot the 
result for the critical concentrations obtained within the momentum-cutoff RG approach, 
Ref. [27]. Both approaches give practically indistinguishable results for the critical fillings 
where antiferromagnetism disappears. Mean-field theory predicts a broader concentration 
range for the stability of antiferromagnetism, see Fig. 9. It was proven by van Dongen 
[47] that for large space dimensionality d ^ 1 and U <^ t the critical hole concentration 
5c — 1 — is reduced in comparison with its mean-field value by a finite factor q^. 
To analyze wether this remains true ior d — 2 we consider the [/-dependence of the ratio 
q{U) = 5c/ S^^. We find that q{U) slightly decreases with decreasing U and it is saturating 
at q{0) ~ 0.4 ±0.025. Surprisingly, the formal application of the results oil/d expansion in 
Ref. [47] to d — 2 gives a close value, g2 0.3. 

We have verified that within the antiferromagnetic phase at t' = the susceptibilities at 
wavevectors Q ^ (tt, tt) are always smaller than the susceptibility xaf, so that the tendency 
towards incommensurate magnetic order is subleading in comparison with commensurate 
(tt, tt) order. The fact that we identify commensurate AF even away from half filling may 
be reconciled with the possibility that the system develops inhomogeneous spin and charge 
structures (e.g. phase separation), as was obtained in the mean-field studies of the Hubbard 
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model [48] and the weak-coupling results in high dimensions [47] . 

Outside the antiferromagnetic region we find the tendency towards dSC as observed 
previously in Refs. [27,34]. The values of the crossover temperature T^g^ rapidly decrease 
away from the antiferromagnetic phase; we show in Fig. 9 the contour lines with ln(i/T^gQ) = 
5, 6, 7. The contour lines with larger T^^q can not be traced within the present RG analysis: 
because of strong fluctuations near the AF phase, the coupling constants reach V^ax = 
ISt before the d-wave susceptibility becomes large. On the other hand, smaller T^^q (and 
correspondingly larger deviations from half filling) are hard to treat, too, because of the 
difficulties with the numerical integrations in Eqs. (15) since the integrands contain sharp 
Fermi functions at small temperatures. 

The half-filled case at different t' was investigated previously within a mean-field analysis 
[49-51], QMC calculations [19,50], and path-integral RG [52]. Different methods predict 
different values of the critical interaction Uc for the onset of antiferromagnetism at fixed 
t'. In particular for t'/t = 0.2 QMC results on an 8 x 8 lattice [19,50] at T = 0.25t yield 
Uc = 2.5t, while path-integral RG [52] gives Uc = 3.4t. The result of the mean-field approach 
for the same t'/t is Uc — 2t [49-51]. We present our phase diagram as obtained from the 
many-patch RG analysis in Fig. 10; symbols show the critical values Uc obtained by other 
methods. As expected and in agreement with previous studies, the critical Uc is larger than 
the mean-field value for all t'. At the same time, the Uc result of the TCRG at t'/t — 0.2 is 
larger than that from QMC calculations, but it agrees well with the path-integral RG result 
in Ref. [52]. 

Again, we find the tendency towards dSC away from the AF region; we show in Fig. 10 
the contour fines which correspond to ln(t/T^g(-,) = 5, 6, 7. Note that from the extrapolation 
of these data to larger U in the paramagnetic phase {U < Uc) we always find \n{t/T^QQ) > 3 
in the weak-coupling regime, i.e. a temperature regime which is far below the accessible 
temperature range in QMC simulations. Therefore, it may be difficult if not impossible 
to observe the corresponding superconducting fluctuations in QMC calculations on flnite 
lattice sizes - at least in the weak-to intermediate coupling regime. 
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C. Ferromagnetic instability 



Now we investigate the ferromagnetic instability, which arises for t' /t close to 1/2. We 
start with t' /t = 1/2, when the dispersion at the bottom of the band at small or ky can 
be expanded as 



i.e. it has extended minima along the lines kx — and ky — (see Fig. 11a) rather than 
a single minimum at the origin, as for t' /t < 1/2. This peculiar flatness of the spectrum 
leads to a square-root divergence of the density of states, p{e) oc at the bottom of 

the band (Fig. lib). Therefore in the low density limit (which is close to a vH band 
filling, since nyn = for t'/t = 1/2), saturated ferromagnetism is expected [23-25]. At 
t'/t — 1/2 the T-matrix approximation [24] predicts rather high critical densities for the 
stability of ferromagnetism, e.g. ric — 0.57 for U — At. For t'/t — 0.45 and U — At the same 
approach predicts ferromagnetism for densities 0.3 < n < 0.5; the smallest value of t' at 
which ferromagnetism can exist was predicted to be {t'/t)c — 0.43 for U — At. The projected 
QMC calculations [25] confirmed the existence of ferromagnetism for {t'/t)c > 0.47. 

The phase diagram obtained within the TCRG approach for t'/t = 1/2 is shown in 
Fig. 12. Similar to the antiferromagnetic instability, mean-field theory overestimates the 
tendency to magnetic order. The result of the T-matrix analysis of Ref. [24] for the critical 
concentration of the stability of ferromagnetism at U = At is marked by a cross. Surprisingly, 
this result is very close to the result of the KG approach. Similar to Ref. [47] one may 
introduce the quantity qpiU) — ric/riiJf^ to measure the deviation from the mean-field result 
at t'/t = 1/2. The analysis of the data shows that qpiU) slightly increases with decreasing 
U and qpiU ^ 0) ~ 0.8. 

We have also explored the possibility for triplet (p-wave) pairing in the vicinity of the 
ferromagnetic phase. Although the p-wave pairing susceptibility is dominant in this region, 
a conclusive low-temperature extrapolation for the inverse susceptibility is not possible. 




(20) 
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Therefore it is not clear whether a finite crossover temperature T*gQ exists. In any case 
the possible values for T*^^ must be significantly smaller than the crossover temperatures 
for d-wave superconductivity. The region where ln(t/Tx) < 8 is shown in Fig. 12, too. 
The growing of the vertices near the ferromagnetic phase results from the triplet p-wave 
superconducting fluctuations, but, unfortunately, the smallness of the temperature crossover 
scale, which is far below the range of applicability of the TCRG method prevents a safe 
conclusion about the possibility of a p-wave superconducting ground state. 

Now we consider the case t'/t < 1/2, which is very different from the above-discussed case 
t'/t = 1/2. The square-root divergence of the density of states is replaced by a logarithmical 
divergence at the energy of the vHS, p{e) oc \n{t/e), while the density of states is finite at 
the lower band edge (see Fig. 13). The phase diagram for t'/t — 0.45 is presented in Fig. 14, 
where we again mark by cross the result of the T-matrix approximation. The ferromagnetic 
region substantially shrinks with decreasing f: it reduces to a narrow density window around 
the vH band filling uvh — 0.465 (the corresponding critical densities are almost symmetrical 
around uvh so that only the region n > nyn is shown). Nevertheless, the ferromagnetic 
region is wider than in the T-matrix approximation. The same tendency is reflected in the 
RG result for the critical value {t'/t)c ~ 0.3 for the disappearance of ferromagnetism at U — 
4t, which is much lower than the result of the T-matrix approximation cited above, {t'/t)c = 
0.43. As well as for the case t' — 1/2 we also find an increasing triplet superconducting 
susceptibility away from the ferromagnetic phase, while the possible corresponding crossover 
temperatures Tpg^ remain undetectably small. 

IV. SUMMARY AND CONCLUSIONS 

We have considered the phase diagrams of the t-t' Hubbard model within two- and 
many-patch RG approaches as shown in Figs. 8, 9, 10, 12, and 14. Instabilities towards 
antiferro- or ferromagnetic order as well as to singlet d-wscve superconductivity are identified 
in different parameter regimes. Near the ferromagnetic region the p-wave superconducting 
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susceptibility is enhanced, but a conclusion about a possible triplet superconducting ground 
state remains elusive. 

At small t' and vH band fillings the antiferromagnetic instability dominates. With in- 
creasing t' antiferromagnetism is replaced by d-wecve superconductivity. At larger t'/ 1 ferro- 
magnetism becomes the leading instability. The tendency towards d-wscve superconductivity 
decreases with increasing U while antiferromagnetism is enhanced. We found that the two- 
patch approach predicts correctly the boundary of the ferromagnetic phase at vH band 
fillings, while it fails to reproduce correctly the boundary between antiferromagnetic and 
superconducting phases at small t' where nearly nesting effects become important. 

Antiferromagnetism at small t' and ferromagnetism at t' /t — 1/2 exist in broad density 
ranges around vH band fillings; the antiferromagnetism remains commensurate in the part 
of the phase diagram where the long-range ordered ground state is expected. The density 
ranges for magnetic order, found from TCRG are substantially narrower (2.5 times for the 
AF instabihty and 1.3 times for the F instabihty at small U) than the corresponding mean- 
field results. 

At half-filling at different t' we find the critical interaction strengths for the antiferromag- 
netic instability. From the present analysis we can not argue, wether the antiferromagnetic 
state we find is metallic or insulating. It was proposed [50] that at nonzero t' there is a finite 
interaction range Uc < U < U'^ for metaUic antiferromagnetism, at C/ > C/^ it is replaced by 
the insulating AF state. On the other hand, the existence of a paramagnetic insulating state 
at larger t' was conjectured in Ref. [52]. Discriminating between these possibilities requires 
the calculation of the conductivity and the Drude weight, for which it is necessary to retain 
the frequency dependence of the vertices. 

The boundary of ferromagnetism at t' = t/2 found from TCRG is surprisingly close to 
the T-matrix approximation result in Ref. [24] at U — 4:t , although the corresponding filling 
is not small and possibly outside the region of the validity of the T-matrix approximation. In 
the vicinity of antiferro- and ferromagnetic phases we found regions with enhanced rf-wave 
and p-wave superconductivity, respectively. Not too close to the antiferromagnetic phase 
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the crossover temperatures for d-wscve superconductivity into the corresponding renormal- 
ized classical regime with exponentially large correlation length can be estimated from the 
extrapolation to low temperatures of RG data for the order-parameter susceptibilities. At 
the same time, triplet p-wave superconductivity in the vicinity of the ferromagnetic phase 
possibly has much smaller crossover temperatures Tpg^ which can not be determined safely 
from the present RG analysis. 

The ferromagnetic phase substantially shrinks for t'/t < 1/2 and the difference to the 
mean-field result increases. For this case, the T-matrix approach underestimates the ten- 
dency towards ferromagnetism. The tendency towards triplet p-wave superconductivity 
in the vicinity of the ferromagnetic phase persists, although its associated temperature 
crossover scale remains very low - significantly lower than for d-wave superconductivity. 

It remains an open issue, how the above results change, when the two-loop corrections to 
the RG equations are taken into account, and how the electronic self-energy evolves in the 
vicinity of magnetic or superconducting instabilities. Another interesting issue for future 
work remains the question whether inhomogenous spin and charge structures are possible 
near half-filling and whether the Pomeranchuk instability may become the leading instability 
for anisotropic extensions of the 2D Hubbard model; if it does it is natural to connect this 
tendency to the stripe pattern formation in rare-earth doped La2-a;Sra;Cu04. Furthermore, 
the tendency towards p-wave superconductivity near the ferromagnetic phase of the t-t' 
model suggests a possible route fo future investigations of the origin of triplet pairing in 
Sr2Ru04. 
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APPENDIX. THE TEMPERATURE CROSSOVER TO THE RENORMALIZED 

CLASSICAL REGIME 



In this Appendix we discuss how the temperature dependence of the magnetic suscepti- 
bihties changes at the crossover to the renormahzed-classical regime. As a first example, we 
consider the 2D ferro- and antiferromagnetic Heisenberg models 

H (21) 

(plus corresponds to the antiferro-, minus to the ferromagnet, J > 0). The susceptibility 
(staggered susceptibihty in the AF case) at high temperatures T ^ J obeys the Curie law 

Xf,af = ^ (22) 

where C — JS{S + l)/3. On the other hand, at temperatures T <g J it was found from the 
RG analysis of the 2D nonlinear sigma model [41,42] 

XF,AF = C'^e (23) 



where 



cxp(27rpjT) AF 
e = Q<; (24) 
{T/jy/^ex:p{2nJS/T) F 

is the correlation length, C^, C are temperature-independent prefactors, and Ps is the zero- 
temperature value of the spin stiffness, which is proportional to the ground-state (sublattice) 
magnetization Sq. Therefore, below the crossover temperature T* ~ 2t^JSq {Sq — S ior a, 
ferromagnet) the susceptibility becomes exponentially large. 

Similar results can be obtained for the Hubbard model within the two-particle self- 
consistent approach [43] . We have 

- X°q(T) .... 

where Q ={n,n) in the AF case, Q = in the F case, and Xq(T) = xo(Q;0,T) is the 
zero-frequency limit of the bare dynamic susceptibility 
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Xo(q, ^^n. T) = Y^ . ^""^^ ■ (26) 



k 

The effective interaction Usp{T) satisfies the self-consistent equation 



U l-C/,p(^)xo(q,^^n,T)• ^ ' 



At high temperatures T ^ t we have UgpiT) ~ ?7 and 

Xf,af^X°q(T)c^^. (28) 

At small temperatures the denominator in Eq. (27) can be expanded for wavevectors q close 
to Q (cf. Ref. [43]) and one obtains 

-o = ^E^-.+ (q_Q)2 (29) 

where Xo = Xo(Q, 0, 0), A = V2xo(q, 0, 0)|q=Q, and al ^ n - n^U,p{0)/2 - V/2 {V is the 
zero-temperature value of the r.h.s. of Eq. (27)), and = [i - [/,p(T)xo(Q,0, T)]/A. 
Momentum integration in Eq. (29) leads to 

e-Qexp{n^a^/2Txo} (30) 

and 

Xf,af XoC^- 

In this case the crossover temperatures Tp^^p ~ tag. Therefore these examples show that 



Xf,af 



C/T T:$>T* 
T<^T* 



and the correlation length ^ oc exp(^T*/T) is exponentially large in the low-temperature, 
renormalized classical regime. 
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Figure captions 

1. The dependence of the van Hove band fiUing on t' /t. 

2. The Fermi surface at van Hove band fiUings: t' — Q and n—1 (sohd hne), t' /t — 0.1 
and n = 0.92 (long-dashed hne), and t'/t = 0.3 and n = 0.72 (short-dashed hne), A 
and B are van Hove points. 

3. The division of momentum space into patches in the two-patch approach. Regions I 
contain the momenta closest to the vH singularity points and k^, regions H contain 
the momenta which are close to the FS but far from vH singularities and regions HI 
contain the momenta far from both, FS and vH singularities. 

4. The four types of vertices considered in the two-patch approach: (a) and (b) correspond 
to exchange and direct scattering between different vH singularities respectively, (c) 
umklapp scattering, (d) intrapatch scattering. The incoming and outgoing momenta 
with equal spin are connected by solid lines inside the vertices. 

5. a) Scaling behavior of the coupling constants at t'/t = 0.1. The solid line corresponds 
to gi, dashed line to g2, dash-dotted line to g^, dotted line to g^. b) Scaling behavior 
of the susceptibilities at t'/t — 0.1. Solid hne corresponds to the AF, dashed line to 
the dSC, and dotted line to the F susceptibility. The interaction strength is U — 2t. 

6. The same as in Fig. 5 for t'/t — 0.45. 

7. Diagrammatic representation for the many-patch RG equations, Eq. (15). Lines drawn 
through the vertices show the direction of spin conservation. Diagrams are drawn in 
the same order as the respective terms in Eq. (15). The cutting dash at the propagator 
lines means the derivative with respect to T (for brevity we indicate only the derivative 
of one of the propagators, the same diagrams with derivatives of another propagator 
are included as well). 
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8. Phase diagram at vH band fillings as obtained from two- and many-patch RG analyses. 
Solid lines correspond to the phase boundaries obtained within the two-patch RG 
analysis. The symbols show the results of the many-patch RG approach: closed circles 
correspond to AF, open circles to incommensurate (tt, tt — 5) order, diamonds to dSC, 
and triangles to the F phase. Long-dashed line is the boundary of ferromagnetic phase 
obtained in Ref. [31] (see text). 

9. Phase diagram for t' = 0. The dashed line is the mean-field phase boundary between 
antiferromagnetic and paramagnetic phases, the solid line is the boundary of the anti- 
ferromagnetic phase obtained from the temperature-cutoff many-patch RG approach. 
The corresponding result of Ref. [27] for the boundary of the antiferromagnetic phase 
is shown by dotted line. The dot-dashed, dot-dot-dashed, and dot-dot-dot-dashed 
lines are contour lines for the d-wave superconductiving crossover temperature into 
renormalized classical regime T^gc = e^^t,e^^t, and e~^t respectively (see text). The 
inset shows the phase diagram in fi-U coordinates. 

10. Phase diagram at n = 1 (half-filled case). The notations are the same as in Fig. 9. 
The cross corresponds to the critical Uc for the stabihty of the antiferromagnetic phase 
at t'/t = 0.2 as obtained from QMC calculations [19], the star marks the result of the 
path- integral RG approach for Uc [52]. 

11. The electronic dispersion (a) and non-interacting density of states (b) at t'/t = 1/2. 

12. The phase diagram for t'/t = 1/2. The long-dashed line MF(F) is the mean- field phase 
boundary between ferromagnetic and paramagnetic phases, the dot-dashed line MF(S) 
is the mean-field result for the boundary of saturated ferromagnetism. The solid line is 
the boundary of ferromagnetic phase obtained from the temperature-cutoff many-patch 
RG approach, the short-dashed line is the contour line above which the maximal vertex 
reach T4iax = 18t at the temperature Tx > e~H. The T-matrix phase boundary for 
the ferromagnetic phase [24] at U — At is marked by cross. The inset shows the phase 
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diagram in //-[/ coordinates., pSC marks the region where the triplet superconduct 
susceptibihty is dominating. 

13. The non-interacting density of states for t' — 0.45i. 

14. The phase diagram for t' /t — 0.45. The notations are the same as in Fig. 12. 
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